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Abstract: Gold nanoclusters have the tunable optical absorption property, and are 
promising for cancer cell imaging, photothermal therapy and radiotherapy. First-principle 
is a very powerful tool for design of novel materials. In the present work, structural 
properties, band gap engineering and tunable optical properties of Ag-doped gold clusters 
have been calculated using density functional theory. The electronic structure of a stable 
AU20 cluster can be modulated by incorporating Ag, and the H0M0-LUM0 gap of 
Au2o-«Ag„ clusters is modulated due to the incorporation of Ag electronic states in the 
HOMO and LUMO. Furthermore, the results of the imaginary part of the dielectric 
function indicate that the optical transition of gold clusters is concentration-dependent and 
the optical transition between HOMO and LUMO shifts to the low energy range as the Ag 
atom increases. These calculated results are helpful for the design of gold cluster-based 
biomaterials, and will be of interest in the fields of radiation medicine, biophysics 
and nano science. 
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1. Introduction 

Gold nanostructures have attracted considerable attention owing to their unique electronic and 
optical properties, as well as their great potential for medical applications [1-3]. Today, gold 
nanostructures are also promising for their possible applications in photothermal therapy, radiotherapy 
and cancer cell imaging because of their good biocompatibility [4-6]. However, a disadvantage of 
photothermal therapy is that the optical absorption of gold nanoparticles, well known surface plasmon 
resonance (SPR), is always located in the visible light region, which limits the wider application of 
photothermal therapy in vivo. To utilize the high optical transmission of biotissue, it is necessary to 
enhance the near-infrared (NIR) optical absorption. In recent years, several groups have developed the 
gold nanorod, nanocage and core/shell structure, which are recognized as the promising materials for 
photothermal therapy [1,7,8]. However, the increasing size would obviously reduce the permeability of 
tumor tissue, which is disadvantageous to photothermal therapy. 

Small gold clusters of typically 1-2 nm are proposed as alternative materials. Sun et al. have shown 
that the Au-SiC>2 clusters can induce the unusual optical transition, and the optical absorption can be 
modulated to the NIR. This is also confirmed by the recent photothermal therapy and drug delivery 
experiment [9]. However, the gold clusters are structurally unstable and have ambiguous electronic 
structure. The procedure of designing gold clusters with NIR optical absorption is influenced 
obviously, because of the long time controversy, by the structure of gold clusters. Recently, the 
obvious progressing has been achieved by density functional theory (DFT) theory calculations and 
optical spectra experiments [10,11]. The common view is that gold clusters prefer the two dimension 
structure in the range of AU4-AU13, while the AU14-AU20 show the cage-like three dimension 
structure [12]. The more complex cage-like AU32-AU38 clusters have been predicted, although the 
experimental result may always be contradicted. The optical absorption of AU2-AU13, AU19, and AU20 
have been calculated by time-dependent density functional theory (TDDFT), while the optical 
transition of AU32 is also focused due to its more stable structure. Notably, the tetrahedral AU20 cluster 
shows a band gap of -1.818 eV, and shows slight NIR absorption [13,14]. 

Doping Au clusters by other metals provides an available route to modulate electronic and optical 
properties [15]. It has been demonstrated that metal atoms, such as Zn and Cu, used to dope Au 
clusters can modify both their structural stability and optical properties due to the delocalizing s and d 
electronic configurations [16]. Bonacic-Koutecky et al. studied the structure and electronic properties 
of bimetallic Ag^Au^ (3 < (M + N) < 5) clusters with DFT calculations [17]. For these clusters, the 
charge transferred from Ag to Au plays a dominant role in the structure of the bimetallic clusters. The 
Au-Ag bond is preferred to the Au-Au or Ag-Ag bond [18]. In addition, the metal-doped Au clusters 
may modulate the HOMO-LUMO gap [19,20]. Thus, we are interested in whether AuAg cluster has 
tunable optical absorption properties. 

Here, we studied electronic structure and optical properties of Ag-doped Au 2 o clusters. The paper is 
organized as follows. Section 2 presents and discusses the results of our calculations. First, we 
investigated the structural properties by analyzing the binding energy. Then we calculated the 
electronic structures, because the optical properties depend on both the interband and intraband 
transitions, which are determined by electronic states. Finally, we analyzed the optical transition in 
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different configurations. Section 3 describes the basic ingredients and details of computational 
methods we applied. Section 4 concludes and summarizes our findings. 

2. Results and Discussion 

2.1. Structural Properties of Ag-Doped Gold Clusters 

Figure 1 gives the calculated ground state geometries of Au20-»Ag„. The stability of gold 
nanostructure is very important for its electronic and optical properties. Thus, the binding energy per 
atom E h (n) of gold clusters is calculated to evaluate the stability of gold clusters, which are defined by 
the following formula: 

ff(Ag) + nE(Au) - rc£(Au„ Ag) 

E b (n) = (1) 

n + 1 

where E(Ag), E(Au), and .E(Au ;i Ag) represent the total energies of the most stable Ag, Au, and AuAg 
clusters, respectively. It is worth pointing out that all of the clusters are found to prefer the lowest spin 
state. We find that the binding energy of AU20 is 2.40 eV, which is very close to the previous 
investigation of gold clusters [21]. It has been demonstrated that the binding energy increases with an 
increase in the size of gold clusters [11]. Increasing binding energy means increasing stability due to 
enhanced core electron configurations. After Ag doping, the binding energy of Au^Agi, AuigAg2, 
Aui 7 Ag 3 , and Aui6Ag 4 are 2.68, 2.68, 2.68, and 2.69 eV, respectively. The increasing doping atom 
induces a tiny effect on the binding energy. It is worth noting that the binding energy of AuAg alloy is 
higher than the Au2o- It shows that Ag atom incorporation can enhance the structural stability. Indeed, 
the Au-Ag bond is stronger than the Au-Au bond and gives an extra o-bonding interaction by the 
overlap between the vacant Ag 4p and valence Au 6s (5d) orbital, which is very similar to the previous 
Ag-doped gold clusters [20]. 

Figure 1. Calculated ground state geometries of Au 2 o-«Ag„. 
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2.2. Electronic Structure of Ag-Doped Aujo Clusters 

Figure 2 shows the density of states (DOS) to reveal the electronic structure of Au2o-»Ag„ clusters. 
The AU20 cluster shows the large HOMO-LUMO gap, which is in good agreement with the other 
computational results [10,21]. The exact band gap of AU20 is 1.47 eV, which is less than the 
experimental data of 1.78 eV (or 1.818 eV) due to the underestimation of electronic states by DFT [22]. 
Meanwhile, the Au d states are dominated in HOMO, and are located in the range of -6 and 0 eV. 
HOMO consists of Au s and d states; the Au d states are dominant. It is clearly seen that when 
n changes from 1 to 4, the DOS also changes. In general, the band gap of Au2o- n Ag„ clusters is less 
than that of the pure AU20 cluster expect for Aui6Ag4. It confirms that the Ag incorporation into AU20 
can induce the obvious effect on gap, which is consistent with the previous results [19,23,24]. The 
exact band gaps of Au2o-nAg„ are 1.42, 1.34, 1.40, and 1.68 eV, which are corresponding to the 
different n values from 1 to 4. The variation of gap can be understood by electronic states. It can be 
seen that the LUMO of Au^Agi, AuigAg2, and Au i7 Ag3 is shifted to the low energy range compared 
with AU20, which can be clearly seen in Figure 2. This shift can lead to the obvious variation of 
electronic properties. As it is incorporated more and more, the p states in LUMO become stronger, 
which leads to the LUMO shift and band gap narrowing. However, the LUMO of Aui6Ag4 shift to high 
energy range, and thus induce the increase of the gap, which should be related to improving structural 
stability and enclosing electronic configurations [25]. Meanwhile, the p and s states are obviously 
enhanced in the LUMO, and these electronic states can have effects on optical properties. 

Figure 2. The partial DOS of (a) Au2o- n Ag„ clusters. 
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2.3. Tunable Optical Properties of Gold Clusters 

For investigating the optical transition of Au2o-«Ag„ clusters, it is necessary to investigate the 
imaginary part of the dielectric function, because it is very important to the optical properties of any 
materials. Regarding AU20 in Figure. 3, there are two main peaks in 82(0)), at 1.79 and 2.51 eV, 
respectively (namely E\, E%), which are very close to the previous results of 1.86 and 2.78 [21]. It can 
be seen in the inset of Figure 3 that direct transition can be caused between HOMO and LUMO. 
Therefore, it can be expected that E\ (1.79 eV) should mainly be caused by optical transitions between 
Au 6s 55% (and Au d, 45%) states in HOMO and Au 6s states in LUMO, which are close to the other 
first-principles evaluation (1.86 eV) [21]. Moreover, it is not far from the HOMO-LUMO gap of 
1.48 eV. In the DOS of AU20, the Au d states produce two peaks, 1.41 and 1.97 eV, which can induce 
some electronic states in the energy level range of -4-0 eV. Thus, the optical transitions of E% can be 
due to the optical transitions between HOMO-1 consist of Au d states and LUMO consist of Au 6s 
(and Au p) states. 

Figure 3. The imaginary part of dielectric function 62(0)) of AU20 clusters. 
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Figure 4 shows the imaginary part of dielectric function S2(co) of Au2o-«Ag„ clusters. Ag 
incorporation induces some obvious variations of optical transitions. Firstly, E\ has gradually 
disappeared, which can be related to the red-shift of £2 and further inhibition of the intrinsic optical 
transition of E\. Secondly, £2 shows the tunable optical properties with the increasing Ag 
incorporation. The E 2 of Au 2 o, AuigAgi, Aui 8 Ag 2 , Aui 7 Ag 3 , and Aui 6 Ag 4 is 2.51, 2.35, 2.25, 2.07, and 
2.01 eV, respectively. To understand these optical phenomena in detail, it is necessary to analyze the 
optical transition by electronic states. 
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Figure 4. The tunable imaginary part of dielectric function S2(co) of Au2o-«Ag« clusters. 
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Figure 5 shows the transition energy level of Au2o-«Ag„ from band structure calculations. The 
optical transition between HOMO-LUMO has been affected by Ag incorporation. The variation in 
energy level can be described in two stages. In the first stage of Au^Agi and AuisAg2, Ag electronic 
states contribute to both HOMO and LUMO. Furthermore, LUMO has slightly shifted to the low 
energy range, which induces the decrease of transition level and can be responsible for the red-shifts of 
E2. In the second stage of AunAg3 and Aui6Ag4, the increasing Ag atom induces more dispersive in 
d states, and LUMO has shifted to the high energy range, which may cause gap variation. It has been 
supposed that the gap widening of Au^Ag4 is important evidence for the enhancement of binding 
energy and structural stability [25]. The heavy doping can have obvious effects on electronic states, 
and further investigation is still interesting. More valuable information is still necessary to probe by 
optical absorption. 

Figure 5. The outline of optical transition of Au2o-«Ag„ clusters calculated by energy level. 
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Figure 6 presents the optical absorption of Au2o-«Ag„ clusters. The dominant absorption around 
450-550 nm can be observed in all Au2o-»Ag n clusters. For AU20, the absorption band of 450 nm could 
be due to the optical transition of Ez, while the 707 nm absorption band is related to the intrinsic 
optical transition of E\. The Ag incorporation induces the red-shift of absorption band (Ez) from 
478 nm (Au^Agi) to 543 nm (Aui6Ag4). We need to consider two possible effects on the optical 
absorption of these Au2o-,Ag„ clusters. On the one hand, the HOMO-LUMO gap can be 
underestimated by the DFT. According to the available optical data of AU20 from the experiment, the 
gap of tetrahedral AU20 is about 1.7—1.8 eV, which is larger 0.23-0.33 eV than the calculated gap of 
1.47 eV. Thus, the actual gap of Au2o-«Ag„ is wider than that of present calculated results. On the other 
hand, the structure stability of Au2o-«Ag„ should also be taken into account. In the geometry 
calculations, the binding energy of Aui6Ag4 is higher than that of the pure AU20 cluster and the other 
AuAg clusters, which indicates the doping feasibility. Therefore, more Ag atom incorporation may be 
promising for further fabrication, NIR absorption and related applications. Zorriasatein et al. 
have proposed that Cu incorporation into Au clusters can modulate the band gap effectively, and the 
results showed that the introduction of Cu enhanced the binding energy per atom compared to 
Au clusters [25]. AuAg clusters also showed good structural stability and optical properties [23,26,27]. 
Thus, it can be expected that Ag incorporation is an effective strategy for modulating optical 
properties of Au clusters. Our results also clearly show that Ag incorporation can modify structural 
stability and modulate the optical properties of AU20 clusters. These methods also have potential 
applications in understanding the optical properties of metal nanoclusters and designing materials for 
photothermal therapy. 



Figure 6. Optical absorption of Au2o-«Ag„ clusters. 
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3. Computational Section 



The calculations are based on density functional theory (DFT) using a plane-wave pseudopotential 
method [28]. We use the generalized gradient approximation (GGA) in the scheme of 
Perdew-Burke-Ernzerhof (PBE) to describe the exchange-correlation functional [29]. Norm-conserving 
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pseudopotential is used to describe the electron-ion interaction. In this code, the plane wave functions 
of valence electrons are expanded in a plane wave basis set, and the use of norm-conserving 
pseudopotential allows a plane wave energy cutoff. Only plane waves with smaller kinetic energies are 
used in the expansion. Reciprocal-space integration over the Brillouin zone is approximated through 
careful sampling at a finite number of k-points using a Monkhorst-Pack mesh [30]. In the present 
system, full electron calculation for the Au2o-«Ag„ clusters is computationally rather expensive, so it is 
better to introduce effective core potentials for the Au and Ag atoms to describe the inner-core 
electrons. Under this approximation, the 5d lQ 6s l outermost valence electrons of the Au atom and 

10 2 

3d 4s outermost valence electrons of the Ag atom are described. It is well known that the interaction 
of a photon with the electrons in the system can be described in terms of time-dependent perturbations 
of the ground-state electronic states. Optical transitions between occupied and unoccupied states are 
caused by the electric field of the photon. The spectra from the excited states can be described as a 
joint density of states between the valence and conduction band. The momentum matrix elements, 
which are used to calculate the 82(0)), are calculated between occupied and unoccupied states which are 
given by the eigen vectors obtained as solution of the corresponding Schrodinger equation. Evaluating 
these matrix elements, one uses the corresponding eigen functions of each of the occupied and 
unoccupied states [22,28,31]: 



where Q is the volume of the elementary cell, v and c represent the valence and conduction bands, 
respectively, k represents the k point, a> is the frequency of the incident light, and u is the vector 
defining the polarization of the electric field of the incident light, which is averaged over all spatial 
directions in the polycrystalline case. The real part of dielectric function E\(co) can be evaluated from 
the imaginary part 62(a)) by the famous Kramer- Kronig relationship. 

o 

The calculations are performed at 20 x 20 A supercell, which contains 20 neutral Au atoms. These 
models of gold clusters refer to the previous work, which has shown that the best structural stability 
and average distance between Au-Au bonding is about 2.73 A [21]. The substitutional method has 
been taken into account in this paper, and Ag atoms are used to substitute Au atoms in vertex site. In 
this way, the Au2o-«Ag„ clusters are built. In order to weaken the interaction with Ag-Ag, the 
separation distance should be as far as possible. We choose the energy cutoff to be 720 eV, and the 
Brillouin-zone sampling mesh parameters for the k-point set are 2 x 2 x 2 for Au2o- M Ag„ clusters. The 
charge densities are converged to 2 x 10~ 6 eV/atom in the self-consistent calculation. And then, the 
Au2o-«Ag„ clusters are optimized with lattice constants and the positions of substitutional atoms. In the 
optimization process, the energy change, maximum force, maximum stress and maximum 

c o o 

displacement tolerances are set as 2 x 10 eV/atom, 0.05 eV/A, 0.1 GPa, and 0.002 A, respectively. 
4. Conclusions 

In summary, a first-principles study has been performed to evaluate the electronic and optical 
properties of Au2o-nAg n clusters in different configurations. The Au2o-jiAg„ clusters show better 
binding energy and structural stability than the AU20 cluster. The increasing Ag concentration can 




(2) 
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induce the HOMO-LUMO gap variation. Subsequently, the optical transition between HOMO-LUMO 
has shifted to the low energy range with the increasing Ag concentration. Tunable optical transition 
has been observed, and was shown to decrease from 2.51 to 2.01 eV with the increase of Ag atoms. 
Our results clearly show that Ag incorporation can modulate the optical properties of AU20 clusters. 
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